Transverse optical plasmons in layered superconductors. 
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We discuss the possible existance of transverse optical plasma modes in superlattices consisting of Josephson 
coupled superconducting layers. These modes appear as resonances in the current-current correlation function, 
as opposed to the usual plasmons which are poles in the density-density channel. We consider both bilayer 
superlattices, and single layer lattices with a spread of interlayer Josephson couplings. We show that our 
model is in quantitative agreement with the recent experimental observation by a number of groups of a peak 
at the Josephson plasma frequency in the optical conductivity of Lai.85Sro.i5Cu04. 



1. INTRODUCTION 

In recent years interlayer Josephson plasmons 
have received considerable attention both theoreti- 
cally and experimentally. So far this attention has 
been restricted to systems consisting of a stack of 
identical 2-dimensional superconducting layers cou- 
pled by identical insulating barriers. In this paper 
we consider the more general problem of ordered and 
disordered superlattice structures. Regardless of the 
underlying microscopic model, in the superconduct- 
ing state the low energy scale properties of a multi- 
layer of Josephson coupled 2-dimensional supercon- 
ducting layers is described by the Lawrence-Doniach 
model[l]. This model has been applied to the case 
of an ordered array of single superconducting lay- 
ers, by Bulaevskii et al.[2], who calculated the longi- 
tudinal dielectric function, and by Tachiki et aL[3], 
who studied the transverse dielectric function. Mi- 
croscopically the interlayer Josephson resonances are 
a result of broken gauge symmetry in a supercon- 
ductor with long range Coulomb interactions [4], and 
belong to a larger class of several types of collective 
modes in layered unconventional superconductors [5]. 
A low frequency Josephson plasmon[6] arises as a di- 
rect consequence of the electrodynamical properties 
of a BCS superconductor in the dirty hmit[7, 8, 9]. 
A different scenario has been discussed by Ander- 
son, in which these modes occur as a direct conse- 
quence of the interlayer pairing mechanism for su- 
perconductivity in the high Tc materials [10, 11]. In 
this spirit we attempt to set the stage for further 



experimental studies of layered superconductors, in 
particular materials with more than one supercon- 
ducting layer per unit cell {e.g. Bi2Sr2CaCu208) 
and materials where the interlayer Josephson param- 
eters are randomly distributed due to disorder {e.g. 
La2-,Sr,Cu04[6, 7]). 

2. RESULTS 

To simplify the discussion we will only consider 
the charge dynamics for fc — > 0, with the electric field 
vector perpendicular to the Cu02 planes, which is 
the relevant limit for a discussion of the c-axis opti- 
cal properties. A discussion of the dielectric tensor 
for general values of k was given by Bulaevskii for 
the superconducting state [2] and by one of us for the 
normal state[12]. Let us consider the relevant Hamil- 
tonian of a stack of superconducting planes coupled 
through the Josephson effect, each with a total area 
A, and with total thickness D of the crystal along 
the z direction. 

Each plane is characterized by a charge per unit area 
Pm and a phase 4'm, where m is the layer index. As 
we are only considering electrical fields perpendicu- 
lar to the planes, the long-range Coulomb potential 
between each set of planes increases linearly with dis- 
tance. The charge dynamics enters via the Joseph- 
son coupling Jm between each set of nearest neighbor 
planes. The equations of motion for a general distri- 
bution of interlayer couplings, and general values of 
wavevector k were derived by Bulaevskii [2]. In the 
present paper we only consider fc — > 0, relevant for 



optical experiments. The solution for general val- 
ues of k will be presented clscwhere[13]. The limit 
k ^ loads to a considerable simplification of tlw^ 
equations of motion, which in this case correspond 
exactly to an equivalent network model where each 
set of two neighbouring layers can be represented by 
a capacitance shunted by an induction. The total 
complex impedance Z of the crystal is the scries resis- 
tance of all local impedances Zm = 4:mdm/{Auiem), 
where we introduce the quantity = 1 — Wj^/w^. 
Here the screened Josephson plasma frequency is 
ATT^^^^e^dmJm, and dm is the 



defined as 

distance between neighboring planes. The expres- 
sion for the macroscopic dielectric function e follows 
immediately after identifying ZA/D with the com- 
plex conductivity 1/ct, and by using e = Ama/w. 
This model can be easily extended to take into ac- 
count quasi-particle currents by shunting each Z^, 
with a resistor in series with an induction, result- 
ing in an expression for e analogous to the Gorter- 
Casimir two-fluid model, with a normal-fluid com- 
ponent an — (Jo/ {I — iijjt). Polarization of the ions 
and transverse optical phonons can be incorporated 
by adding separate oscillators to the expression for 
e. The expression for e becomes 



pole at (jJt, which corresponds to a transverse opti- 
cal plasma oscillation, with oscillator-strength St = 
€00(1 — u)j/u)^){l — w^/wy). Photons with frequen- 
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where Wm = dm/D is a weighting factor proportional 
to the distance between the layers. The total dielec- 
tric function is a sum over response functions of sets 
of two coupled planes for which we obtain 
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3. BILAYER SUPERCONDUCTORS 

Let us now consider a crystal which has two su- 
perconducting layers per unit cell. We denote the 
Josephson coupling constants as I and K, the dis- 
tances between the planes as dj and dK^ and the 
screened Josephson plasma frequencies as uji and ujk- 
Using Eq. (3), and approximating the effect of ionic 
and vibrational screening with a constant Cqo, the 
dielectric function is 



(3) 



where we introduced lo^ = {djuj^ + dK<jj'j)/{di + 
dx). Interestingly, there are now two longitudi- 
nal plasmons. Also the dielectric function has a 



Figure 1: Reflectivity, Ree, optical conductivity and 
loss function calculated for a hypothetical bi-laycr su- 
perconductor, assuming wj = 5cm~^,a;T = 26cm~^, 
and wk = 30cm~^. The chained curves were cal- 
culated assuming no normal-fluid component in the 
conduction, the solid curves were calculated assum- 
ing a parallel conductance of 2 S/cm. 

cies within the Reststrahlenband lot < lu < lok, and 
with E polarized along z can not propagate inside the 
crystal and are totally reflected at the interface. Such 
a reflectivity spectrum of the crystal-face perpendic- 
ular to the superconducting layers is displayed in Fig. 
1, along with the dispersion relation. We see, that 
a lifting of degeneracy between the two Josephson 
frequencies and ujk should be easily recognizable 
in the reflectivity curve, even if the features in a are 
blurred due to parallel conduction of quasi-particles. 
The compound Bi2Sr2CaCu208 may exhibit such be- 
haviour. There a resonance found at 5 cm~^[14, 15] 
is probably the lowest of the two Josephson plas- 
mons. The higher plasma frequency has not yet been 
observed with infrared spectroscopy down to w 30 
cm~^[16] and can be situated anywhere between 5 
and 30 cm~^. 



4. SINGLE LAYER SUPERCONDUCTORS 
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Figure 2: Experimental loss function (upper panel) 
and conductivity (lower panel) of Lai.85Sro.i5Cu04 
at 10 K (open circles), 15 K (close circles), 20 K 
(open squares), and 27 K (closed squares). The solid 
lines are a fit with the same relative gaussian distri- 
bution with a FWHM which is 20% of the centerfre- 
quency. 

We now apply Eq. |l| to a stack of one supercon- 
ducting layer per unit cell, but with a distribution 
of the interlayer coupling around a center value. As 
was the case for the ordered superlattices discussed 
above, the material is effectively homogeneous vis 
a vis the optical response of a lightbeam reflecting 
on the ac-face of the crystal, as long as the length 
scale of variations of is much smaller than the 
wavelength of the electromagnetic radiation. Let 
us now consider the case where a spread exists of 
the Josephson frequencies around a center frequency. 
The source of such could be the local variations in the 



potential barrier, e.g. due to random substition of 
La'^^ with Sr^+ ions in the tunneling barrier between 
the superconducting Cu02 layers in La2-xSra;Cu04. 
The expression for e then becomes 
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where p{ujj) is the normalized distribution function 
of the screened Josephson plasma frequencies. The 
distribution p{ujj) should decay sufficiently rapidly 
for ujj 0. Physically this corresponds to the con- 
dition that the weakest Josephson junctions in the 
stack of layers should still have a finite critical cur- 
rent (or ujj value). 

A peak in the optical conductivity which coincides 
with the longitudinal Josephson plasma has been re- 
ported by at least four different groups[6, 17, 7, 18], 
but has to our knowledge not been explained in a 
satisfactory way. In Fig. 2 we display both the ex- 
perimental loss function and the conductivity for the 
region near the Josephson plasma frequency. These 
data were adopted from the paper by Kim et al. 
So far we had the data of Gerrits et al. and Kim 
et al. at our disposal, and checked the validity of 
our approach for all temperatures reported for these 
two sets of data. The agreement with the formalism 
outlined above is very satisfactory at all tempera- 
tures. To demonstrate this result we display in the 
top panel of Fig. 2 the loss function at 10 K of Kim 
et al.. In the lower panel the optical conductivity is 
displayed. Clearly both the loss function and cr peak 
at the same frequency of 48 cm^ ^ . If wc assume that 
the 'intrinsic' value of the c-axis conductivity is 1.8 
S/cm in this case (dashed line in the lower panel), 
and assuming that there is a single value of wj of 
48 cm~^, we obtain the dashed curve for the loss 
function in the top panel. The finite width in this 
case is due to the presence of a finite background 
conductivity, and is given by 2/r = Sna-n/coo. With 
(T„ = 1.8 S/cm and e^c ~ 22 we obtain c/t = 27r • 10 
cm^^. In other words: Experimentally the broaden- 
ing of the plasmon-peak can not be explained as the 
usual lifetime effect. If we now assume that in addi- 
tion ujj has a Gaussian distribution, we obtain the 
solid and chained curves in Fig. 2 for the loss func- 
tion and the conductivity using Eq. ^ with a FWHM 
of 11 and 17 cm~^ respectively. Clearly the choice 
of FWHM which gives good agreement for the loss 
function, also explains the anomalous peak in the 
conductivity both qualitatively and quantitatively. 



5. CONCLUSIONS 

We discussed the possible existance of transverse 
optical Josephson plasmons in layered superconduc- 
tors with 2 or more layers per unit cell. A related 
-but not similar- effect effect may have already been 
observed in La2-2:Sr2:Cu04 due to variations in the 
interlayer Josephson frequency, possibly due to the 
random substitution of La with Sr. It should be pos- 
sible to design and grow artificial superconducting 
multilayers with optical properties taylored to spe- 
cific needs in the THz range. Because they provide 
an effective channel for coupling electromagnetic ra- 
diation to the current perpendicular to the super- 
conducting layers, transverse optical Josephson plas- 
mons may become useful for the detection and emis- 
sion of far-infrared radiation in the THz range. 
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